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In this paper we study a methodical problem related to the magnetic scenario recently suggested 
and initiated by the authors [1] to understand the strongly coupled quark-gluon plasma (sQGP): 
the electric flux tube in monopole plasma. A macroscopic approach, interpolating between Bose 
condensed (dual superconductor) and classical gas medium is developed first. Then we work out a 
microscopic approach based on detailed quantum mechanical calculation of the monopole scattering 
on electric flux tube, evaluating induced currents for all partial waves. As expected, the flux tube 
looses its stability when particles can penetrate it: we make this condition precise by calculating the 
critical value for the product of the flux tube size times the particle momentum, above which the flux 
tube dissolves. Lattice static potentials indicate that flux tubes seem to dissolve at T > Td iaso i ution m 
1.3T C . Using our criterion one gets an estimate of the magnetic density n w 4.4 ~ 6.6/m -3 at this 
temperature. 



I. INTRODUCTION 
A. Electric-magnetic duality and monopoles 

Studies of magnetically charged objects resurface in 
theoretical literature regularly since Maxwell's time, as 
soon as some new developments call for their application 
in different fields. Three of these are especially impor- 
tant: (i) the celebrated Dirac quantization condition [2], 
(ii) the idea put forward by 't Hooft and Mandelstam 
[3] to view the QCD confining vacuum as a "dual su- 
perconductor", and (iii) Seiberg and Witten [4] studies 
of ftf = 2 SUSY gauge theories, identifying magnetically 
charged degrees of freedom and the usefulness as well as 
necessity of jumping from "electric" to "magnetic" lan- 
guage whenever appropriate. There are excellent reviews 
covering historic development of these ideas, see [5-7], 
and (for more recent progress) [8,9], and also the excel- 
lent book [10] by Shnir. 

In the context of QCD phase diagram (finite tem- 
perature T - baryonic chemical potential /i) the issue 
of electric-magnetic duality were discussed in our pre- 
vious paper [1] (to be referred below as LSI). Brief 
summary can be made by noticing that since electric 
and magnetic objects generally repel each other, their 
simultaneous presence results in dominance of one over 
the other. In particular: (i) at high-T the Quark-Gluon 
Plasma (QGP) is dominated by electric objects (gluons 
and quarks); (ii) at high-/i and low T QCD is in a color 
superconductor phase, dominated by electrically charged 
diquarks and confining magnetic objects; while (iii) the 
vacuum is presumably a dual (magnetic) superconductor 
confining electric objects. 

The most important point made in LSI is an intro- 
duction of new "electric-magnetic equilibrium line" , on 
which both couplings (electric and magnetic), densities, 
screening lengths etc are about equal, and thus the vol- 



ume is somehow shared equally. The Dirac condition 
then demands that at the equilibrium line both electric 
and magnetic "fine structure constants" 1 are 



Furthermore, as pointed out by 't Hooft, in nearly all phe- 
nomena the relevant effective coupling is not e, g them- 
selves but A e = e 2 N c or X g — g 2 N c , further increasing 
with the number of colors N c . Thus one naturally is lead 
to the conclusion that the QCD plasma at e/m equilib- 
rium conditions must be strongly coupled. In LSI this 
was proposed to be a possible explanation of why exper- 
iments at RHIC have indeed found a strongly coupled 
QGP, see [11-14] for discussion of data and other theo- 
retical ideas. 

The QCD vacuum and deconfining transition was at- 
tempted to be described by an ensemble of self-dual 
dyons in [15] [16]. The presence of magnetic monopoles, 
both below and above the deconfinement phase transition 
was studied in vast literature on lattice gauge theories, 
see e.g. [17]. Our previous paper LSI has also studied in- 
terrelation of electric and magnetic quasiparticles numer- 
ically, but using much simpler classical tool - the Molec- 
ular Dynamics (MD). While lattice study uses Euclidean 
time formalism and thus is restricted to thermodynami- 
cal observables only: we on the other hand were able to 
do real-time simulation and calculate such kinetic prop- 
erties as viscosity and diffusion constant, see [1] for more 
details. 



1 Note in this work we change units from LSI and use cou- 
plings normalized as standardly done in gauge theories, with 
a — g 2 /4ir in the Coulomb law, i.e. Heaviside-Lorentz units. 



B. Electric Flux Tube Formation in Magnetic 
Plasma 

After this brief introduction, let us turn to the subject 
of the present paper, the flux tubes. In the usual electric 
superconductor the corresponding solution of Ginzburg- 
Landau equations was first found by Abrikosov [18], 
which was later revived in field theory [19] and became 
known as Abrikosov-Nielsen-Olesen (ANO) vertex. In 
the "dual superconductor" picture of the QCD vacuum 
properties of the QCD confining string and the resulting 
heavy quark potentials have been discussed extensively: 
see e.g. reviews by M. Baker [20], and more recetnly by 
G. Ripka [21] (with exhaustive list for further references). 

Lattice studies (e.g. [22]) provided substantial support 
to these works. Flux tube behavior at finite T was also 
extensively discussed: in particular Polyakov [23] has 
shown how exponential growth of flux tube entropy leads 
to vanishing of the effective tension in free energy F(r, T) 
and Hagedorn-like phase transition. This scenario would 
predict gradual deconfinemcnt with the string tension 
vanishing at T c : in fact for N c > 2 it jumps to zero. 
Deconfinemcnt transition for various number of colors 
iV c was studies in detail: see e.g. [24] where N c up to 12 
was studied. Working with metastable "overheated" con- 
fined phase it was found that the Hagedorn-like transition 
(at which the string tension of the free energy vanishes 
<r(T) — ► 0) can be approximately located at a universal 
(N c independent) T H /T C = 1.116(9). 

Heating usual superconductors above the critical tem- 
perature destroys not only the condensate but also 
Cooper pairs themselves. Although normal (metallic) 
phase is a plasma of electric objects (electrons), but their 
characteristic momenta p ~ pp are orders of magnitude 
larger than momenta of Cooper pairs, thus there is no 
analog of Abrikosov vortexes in the normal phase. This 
does not happen because presence of a quantum conden- 
sate is that necessary for flux tube's existence: a coun- 
terexample can be provided e.g. by quite spectacular 
magnetic flux tubes in solar plasma 2 . Whether charges 
are Bose-condensed or not, their scattering on a flux tube 
may provide a pressure which may lead to its stabiliza- 
tion. It is just a matter of certain quantitative condition 
for tube stabilization being met. 

The questions to be addressed in this work is whether 
QGP is like an electric plasma in a metal, without mag- 
netic flux tubes, or like other plasmas which have them? 
What exactly are the necessary conditions for a flux tube 
formation? Below we will ignore electric quasiparticles 
which would induce screening/termination of electric flux 



lines [25] and consider purely magnetic plasma. We will 
perform quantum-mechanical study of monopole scatter- 
ing on the tube and examine their back-reaction to the 
tube field through the associated magnetic current. This 
will answer these questions. 

But before we do so, let us explain few important is- 
sues classically, related to the very essence of the electric- 
magnetic competition, i.e. "expulsion" of sub-dominant 
component into flux tubes and their stabilization. A full 
quantum mechanical calculation will be presented in Sec- 
tions III-V. 

We first start with an electric charge e being placed 
within a free gas of monopolcs with mass M and charge 
±g. The monopole gas should be neutral, i.e. with equal 
number of positive and negative charges. We emphasize 
in advance that monopoles with either signs have the 
same effect: this will be seen in the appearance of g 2 
rather than ±g in the final results. 

At a distance R from the charge (see Fig.l left), the 
(unmodified) electric field = 4 ^ R i R will stir the mag- 
netic monopoles into Larmor motion with radius r^. As 
Poincare has shown [5,6,10] a century ago, the radius 
shrinks near the charge, restricting the motion to a cone 
- a small patch of the whole space solid angle. The 
cone angle is determined by (with v t the monopole ve- 
locity transverse to f) 



cot9 



(ge)/4irc 
Mv t r 



(2) 



The numerator is precisely the field angular momentum 
of a charge- monopole pair Lem = (ge)/47rc as first com- 
puted by J. J. Thompson in 1896 [5,6,10], while the de- 
nominator is the monopole's kinetic angular momentum 
L v — Mv t r with respect to the origin. The above for- 
mula, rewritten as cot6> = Lem/L v , reflects the inter- 
play between angular momenta of the electromagnetic 
field and of the particle motion. Though superficially L v 
is defined through v t and r, it is actually a conserved 
quantity uniquely related to the cone angle 9, see [10] for 
detailed discussion. 

In turn, these monopolcs form loops of magnetic cur- 
rent griLL v /Mr (til their density) on the cone. The di- 
rection of the current explains the sign of induced electric 
dipolc 3 . Using dual Maxwell's equation yx£ = — ^ Jm, 
one finds that such electric dipolc is opposite to induced 
dipoles in dielectric, so in this sense it is an anii-screening 
effect. The charge repels such a dipole: thus monopole 
will fly away from the charge. 

To make this statement quantitative, let's calculate the 
curl of magnetic current around R. To do that we need to 
require that the Larmor circle to be fairly small, for two 



2 They have very large fluxes and sizes, and thus a macro- 
scopic theory - magnetohydrodynamics - can be used for their 
description, which unfortunately it is not applicable in our 
case, for microscopically small tubes. 



3 Note that although monopoles with ±g rotate in opposite 
directions, they produce currents of the same sign, so it is not 
necessary to distinguish them. 



important reasons: (i) if it is not small then one has to 
take into account the variation of electric field strength 
which will warp the circle; (ii) a small Larmor radius 
enables one to approximate the V x J M by integrating 
Jm along the circumference and dividing it by the area of 
the circle. Small Larmor radius means small angle 9, 
i.e. L v << Lem- Density of monopoles with angular 
momentum L v at R, is related to total monopole density 
n by til — /lii, with some function of R and 9. In 
such case, the result comes out as: 



(V x J M ) R 



gn L (L v /Mr){2irr L ) 



V(c/\1)E* 



R 



(3) 



Here Xl = (Mc 2 j g 2 n) x l 2 is the London penetration 
length. Interestingly enough one arrives at the second 
London equation with an modification coefficient 



V = 2f L cos9 = 2f L 
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(4) 




FIG. 1. (color online) Schematic demonstration of an- 
ti-screening effect and flux tube formation in magnetic 
medium, see text for detailed discussions. 

To end this discussion, we emphasize again its main 
point: motion of plasma monopoles can function simi- 
larly to condensed monopoles, under the important con- 
dition that the field is strong enough, which specifically 
means that monopoles' bombarding angular momentum 
L v is much less than the field angular momentum Lem- 
The effect discussed here and the well-known Meissner 
effect of (dual) superconductor share the same mech- 
anism, namely magnetic particles are scattered by the 



Lorentz force, inducing currents which anti-screen the 
field. On the other hand, there is a big difference be- 
tween the two. Meissner effect is only present in super- 
conductors, which expel arbitrary weak magnetic field. 
Plasmas do not have it, and thus distributed weak fields 
can be present in the bulk. However, strong enough fields 
can be expelled by the plasma into flux tubes, which are 
in principle metastable. Solar plasma (for example) has 
both weak distributed magnetic fields as well as magnetic 
flux tubes (visible in telescopes as substructure of famous 
solar "dark spots" ) . 

Let us now further simplify the problem, by removing 
the electric charges to infinity and leaving only a flux of 
electric field, confined in a flux tube (see Fig.l right). 
Again, monopoles from outside which move into it are 
turned away by Lorentz force and leave. Although their 
energy remains unchanged, the momentum is changed, 
which means that there is a constant pressure acting from 
the monopole ensemble on the tube. In essence, it is just 
dual to the Meissner effect of electric superconductor, in 
which the magnetic field gets pushed away. 

The situation however is only simple if the strength 
of the field in the tube is such that particles penetrate 
only small part of its radius. (One can then further sim- 
plify the problem into flat surface, with effective current 
floating in wall separating field-free and field regions, as 
is the case for magnetic flux tube in solar plasma.) As 
seen from the Fig.l above, if particles can penetrate into 
the flux tube beyond its center, they start generating a 
counter-rotating current which eventually destroys the 
tubes. To further shed light on how the microscopic flux 
tube may or may not exist, more careful analysis of the 
induced current during scattering of these bombarding 
monopoles will help. 

The lower right picture (transverse projection of the 
upper) shows a few trajectories (with different ending 
points A, B, B' ,C, C") in the constant field E region 
(within tube radius R), starting from the same initial 
point / (at the bottom) with same velocity v (thus curled 
with same Larmor radius = Mvc/gE), yet with differ- 
ent impact parameter b as they are aiming at different an- 
gles. The impact parameter is related to the monopole's 
incoming angular momentum by L v = ±Mvb depending 
on whether the velocity orients toward left or right at the 
initial point (noting the positive z is pointing into the 
page in the figure). For example, the red curve(/ — * A) 
is for 6 = and L v = 0, the blue solid/dashed curves 
is for certain nonzero b and same \L — v\ yet opposite 
signs (with / — > B positive and / — > B' negative), and 
the green solid/dashed curves for even larger b and \L V \ 
(with / — > C positive and I —> C negative) . 

We now consider the currents Jb produced by various 
trajectories. In particular let's examine how the com- 
bined current jb of two trajectories with same b and 
±|£„| changes with b. The important observation is the 
following: trajectories with small b or small \L V \, like 
I — > A, contribute counterclockwise currents, while tra- 
jectories with large b or large \L V \, like / — > C, C", con- 



tribute clockwise currents, and there is a critical b or 
\L V \ (which is precisely the solid blue curve / — > B go- 
ing right through the center of tube) beyond which the 
current inversion happens. Some simple algebra leads to 
the following critical angular momentum of the current 
inversion: 

it i i g®E 



2ttc 



(5) 



with $e — EttR 2 the electric flux. By interpreting the 
right-hand-side as the electric-magnetic field angular mo- 
mentum in this cylindrical setting, we simply have critical 
momentum \L C \ — Lem which coincides with the anal- 
ysis in the previous example. Very importantly, strong 
electric field means large \L C \ and stable flux tube, while 
weak field (with small \L C \) prefers becoming diffusive in 
the bulk rather than expelled into flux tube. 

As we will show later in the paper, this current in- 
version phenomenon is very important. The counter- 
clockwise currents (from small L v ) strengthen the origi- 
nal field 4 , while the clockwise currents (from large L v ) 
weaken them. Thus the current inversion is like a 
"para/dia-electric" inversion, in macroscopic language, 
and it kills the flux tubes. 

To summarize the lesson from this classical example, 
the value of the angular momentum plays essential role 
in the monopole scattering by the flux tube. If particles 
have typical momentum p and the radius of the tube is 
R the angular momentum is L ~ pR. When L is small 
or equivalently the electric field is strong, the motion is 
still basically radial and the pressure argument works. 
In the opposite limit of large L ~ pR » Lem or weak 
field, the induced currents have both signs and cancel 
each other, and there is no reason for flux tube to ex- 
ist. Thus there must be some critical value of pR above 
which there is no flux tube solution, depending on exact 
magnitude of currents induced in channels with different 
angular momentum, to be evaluated quantum mechani- 
cally below. 

The rest of the paper is structured as following: by 
devising and solving a generalized London's equation, 
we first show in Section. II how electric flux tube solu- 
tion could follow from " macroscopic" electrodynamics in 
medium beyond superconductor; in Section. Ill we outline 
our self-consistent treatment of flux tube starting from 
"microscopic" level; the quantum mechanic scattering of 
single monopole in flux tube field will be exactly solved 
with analytic wave functions presented in Section. IV for 
both non-relativistic and relativistic cases; we then pro- 
ceed in Section. V to impose self-consistent condition and 



4 Note again that in the (dual) superconductor case [21], 
the Abrikosov vortex is exactly supported by supercurrent of 
scattered condensate in lowest possible angular momentum, 
namely L v = channel. 



find flux tube size in a thermal medium; application of 
our results to sQGP problem is discussed in Section. VI; 
and finally Section. VII is for conclusions. 



II. ELECTRIC FLUX TUBE: MACROSCOPIC 
APPROACH 

Borrowing wisdom from electrodynamics of a super- 
conductor and being motivated by the "modified Lon- 
don" relation (3) we discussed above, we find in this sec- 
tion solutions to macroscopic electrodynamics equations 
of London's type. 

Our generalized (dual) version of the second London 
equation reads: 



V x Jm = 



A 2- 



-E 



(6) 



Any constant coefficient could be absorbed in a re- 
definition of A. When combined with one of the (dual) 
Maxwell's equations y x E = —^Jm, it yields the equa- 
tion for the electric field 



V 2 £ = 



E 



(7) 



Macroscopic parameter k characterizes how the electric 
field gets modified by the magnetic medium, k = is the 
London limit (appropriate for the medium being a dual 
superconductor in extremely typc-I regime), while k = 1 
corresponds to the classical monopole gas (as discussed 
in preceding section). Intermediate values of k are sug- 
gested as an interpolation between the two limits, say 
to describe a medium having both Bose condensed and 
non-condensed components. 

Our setup corresponds to cylindrical flux tube (see 
Fig.l right), with E = E(r)z in coordinates (r, <j>, z). The 
total electric flux is <&e — E(r)2irrdr. 

The solution for any n < 2 is given by 5 



E(r)= f K % K 



:(t) 



(2-k)/2 



2 — k X J 
f K = l/[(2- K )^r[2/(2- K )] 2 ] 



(8) 



with A'oM, r[a;] being the Bessel and Eulcr Gamma 
functions. If a function F [r/X] is a solution to Lon- 
don eq.(7) with k = 0, then the function F K [r] oc 

- F o[2^( r / A ) (2_K)/2 ] is a s °l ution t0 the modified eq.(7) 
with any n. The normalization constant follows from the 
total flux value. 

At large distance, r — * oo, the electric field 



E ■ 



exp[- (r/A)* 1 -"/ 2 ) /(l - k/2)] j (r/A)^/ 4 



(9) 



3 For k > 2 the boundary condition couldn't be satisfied. 



vanishes quicker than exponential, leaving most of the 
flux within r ~ few A. The smaller is re, the thinner is 
the flux tube. 

For this flux tube solutions the "string tension" - the 
energy per unit length along z is 



E(r)\ , $2 



2 (6-k)/(2-k) . ^ _ k )(2+«)/(2-k) 

r[(6- k)/(4- 2k)] • r[2/(2- re) 



(10) 



7^ is a rapidly decreasing function of re, and in particular 
T K=0 /T K=1 = 3 which means that (for fixed A) a quantum 
condensate expels the electric flux into a flux tube with 
the tension three times larger than a classical monopole 
gas does. 

Let's summarize the physical picture as well as limita- 
tions of the established solution. Any magnetic medium 
generically expels the electric field, as monopoles are 
back-scattered off it, so there is possibility for flux tube 
formation. But different media do this job with distinct 
efficiencies, leading to flux tube (if there is any) with 
different tensions. The electric flux tube solutions are 
rather simple: they describe the problem in terms of two 
macroscopic properties of magnetic medium, namely re 
and A. There is however an important limitation. The 
macroscopic approach is suitable only if the electric field 
strength (or the electric flux here) is large, as the de- 
tailed analysis in the introductory examples has shown: 
we repeat that strong field makes small Larmor radius of 
monopoles, thus the flux tube is a "macroscopic" object 
and scattering of monopoles happens basically on the sur- 
face. If it is not so, the validity of eqn.(6,7) upon which 
the solution is based is no longer justified. The intermedi- 
ate case between diffusive weak field and macroscopically 
strong flux tube requires a microscopic approach, to be 
discussed in the following sections. 



III. THE MICROSCOPIC APPROACH 

Starting here and following in subsequent sections, we 
will pursue a fully quantum mechanical microscopic ap- 
proach. Let us first describe our strategy and approx- 
imations made in this section. The main one is that 
mutual interaction among monopoles will be neglected, 
as it has been argued that magnetic sector of sQGP at 
just above T c is very weakly coupled, see [1] for more de- 
tails. What's more, if one assumes the monopoles are of 
't Hooft-Polyakov type [26] , the Coulomb interaction be- 
tween monopoles may be largely cancelled (and exactly 
cancelled in the BPS limit [27,28] for static monopoles) 
by scalar/Higgs exchange. But the Lorentz force from 



electric field 6 cannot be cancelled and this is the only 
interaction of monopoles relevant to our approach. The 
single monopole scattering on a flux tube will be treated 
quantum mechanically Both non-relativistic and rel- 
ativistic cases will be analyzed: there are evidences 
that monopoles in sQGP are semi-relativistic, e.g. with 
M~2T [29,30]. 

In the following sections we will go through the three 
steps below: i) first assume existing flux tube of certain 
size R, ii) then figure out in great details how individ- 
ual monopole from medium will be scattered off it and 
generate some magnetic current, and iii) finally use the 
dual Maxwell's equation relating the electric field and 
magnetic current to obtain a self-consistent equation de- 
termining the value of R (and thus string tension a) as 
a function of medium parameters T, n, M . Below we ex- 
tend the description of strategy a bit more step by step. 

i) For our purpose the electric flux tube with flux <j># 
and size R is described in cylindrical coordinate (r, </>, z) 
by the following field: 



E 



-{ 



Ejz = <S>e/(kR 2 ) z , r <R 



, r> R 

The corresponding dual vector potential reads: 



C = C A 6 = 



$g r_ 

2ttR R 

_*e_ R 

2ttR r 



r <R 
r>R 



The string tension is given by 

a = $|/(27ri? 2 ) 



(11) 



(12) 



(13) 



According to Dirac quantization, the flux can be normal- 
ized via (g$ e) / {4nTic) — d/2. While the results obtained 
below can be used for general d, we are particularly inter- 
ested in d — 2 as is true for adjoint monopoles in sQGP 7 . 

ii) A monopole moving in such a field is governed by 
the following Hamiltonian: 
in the non-relativistic case 



H NR = (p + 9 -C) 2 /(2M) 

c 



while in relativistic case it is 



Ha = \ P-"c) c 2 + M 2 c* 



(14) 



(15) 



6 Although the monopoles in sQGP are built out of non- 
Abelian fields A^, each type of monopoles only interacts with 
electric field projected into their corresponding U(l) (see e.g. 
[10,21]), so the Maxwellian field description still holds. 

7 There are strong evidences from lattice study of high-T 
magnetic QCD which supports the idea that monopoles in 
QGP have such charges that d = 2 and their total numbers 
scale as N% — 1, see e.g. [31] 



The conserved quantities are 

1) total energy e; 

2) longitudinal momentum p z ; 

3) hence one can use longitudinal energy e z = p 2 /2M and 
transverse energy e t = e — e z to be conserved separately; 

4) the angular momentum L z = r {p^ + - C^). 

The conservation of both e t and L z implies that the 
monopole is rejected back (unless L z = 0) when it ap- 
proaches the center of the tube, due to generic centrifugal 
barrier ~ L 2 /(2Mr 2 ) which dominates at small r. 

iii) The dual Maxwell equation, yx£ = — i J M , in cylin- 
drical setup being 

can be integrated in r 

E{r = 0) - E(r = R) = -- [ jfjdr (17) 

c Jo 

The flux tube may presumably be approximated by a 
constant E inside certain radius r < i?_ and zero out- 
side r > R + , with smooth interpolation in between. As 
an approximation in step(i) we have used step-like elec- 
tric field, neglecting the difference between R±. The ad- 
vantage is that monopole motion in such field can be 
calculated (step (ii)). This shouldn't be a serious issue 
as we expect R + — R- << R. Thus we take E(r — 0) as 
the constant field strength Ej = ^e/(^R 2 ) within tube 
and send E(r = R) to zero, obtaining the equation to be 
used in later section: 

1 r R 

Ej = <S> e /(ttR 2 ) = — / Jffdr (18) 
c Jo 



IV. QUANTUM MECHANICAL MOTION OF A 
SINGLE MONOPOLE 

Quantum mechanical motion of single monopole is de- 
scribed by wave function W which is a scattering solu- 
tion to TL^i = with H from eq.(15). Making use 
of conserved quantities, we may decompose the wave 
function into $ = f(r)e" ncf 'e lKzZ , with energy 8 e = 
e z + e t = (HK Z ) 2 /2M + (hk) 2 /2M and angular momen- 
tum L z = mh. Let's first introduce several parameters 
involved later in the solution, including: 



Here we first deal with non-relativistic case, while in the 
last subsection the treatment will be generalized to relativistic 
case which turns out to be rather straightforward. 



v = m + d (19) 
7 = 1 + H (20) 
a = (kR) 2 /(4d) - m/2 (21) 

We repeat that d — e) / {2nTic) tells how much flux 
is going through the tube. The meaning of v can be 
explained as follows: it is quantized (integer- valued) form 
of a relation between velocity, canonical momentum and 
dual field rav — p + gC projected (via their cross-product 
to r) to angular momenta. Classical path / — > A in Fig.l 
which has velocity at large distances directed to the tube 
center corresponds to v = 0. The m = channel is 
the one corresponding to / — > B path: it goes through 
the center because it experiences no centrifugal barrier 
~ m 2 jr 2 . As we will see below, this correspondence will 
explain the signs of the currents, generated in each partial 
waves. 

The Schrodinger equation can then be reduced to the 
following cylindrical radial equation 

;s(-^) + [*"-"-/] A--" « 22 » 

The effective potential takes the form: 
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FIG. 2. (color online) The effective potential V e ff as a 
function of r/R for v = 0(red), v — ±2(blue solid/dashing), 
and v = ±5 (green solid/dashing). 

The equation could be exactly solved both inside and 
outside the tube, as shown separately below: 
i) r < R (inside), the solution is given by confluent hy- 
pergeometric function i_Fi[x]: 



dr 2 

7/2 -a, 7,-^2 



(24) 



ii) r > i? (outside), the solution is expressed by two Han- 
kel functions H^ ,2 \x] with proper phase shift S k y. 



(ii) 
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Finally the two functions should be connected 
smoothly at r — R, which determines: 
the normalization constants A, B satisfying (with 
J„ [x] , Y v [x] Bessel functions) 



TZab= 
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k.v 



k.v 



B, .. 

v {cos 5 um) J v[kR] - {sh\8 k ,u)Y v [kR] 



e -d/2 d ( 7 -i)/2 lFl [ 7 / 2 - a, 7, d] 
and the phase shift 5k, m being 

' J u+1 [kR] - GJ u [kR] 



8k,v= arctan 
G 



Y v+1 [kR] - GY v [kR] 
v-{-)-l-d)-(l-2ah)-d-F / (kR) 



(27) 



F= 1 F 1 [ 7 /2-a+ 1,7+ l^]/ 7 ifi[7/2- a, 7, <*| 

However exceptions to eq.(26,27) can occur when it so 
happens that ii*i[7/2 — a,-f,d] = 0. In such situation 
the alternative equations are the following: 



T^ab — 



A 



k.v 



(28) 



e* 6 *." (kR) [{ca&6 k , v )J v+1 [kR] - (sm8 k , v )Y v+1 [kR]] 
e -d/2 d ( 7 +i)/2 ( 2a /7 - 1) 1F1 [7/2 - a + 1, 7 + 1, d] 



Sk.u— arctan 



MkR] 
Y u [kR] 



(29) 



The coefficient B k , v should be determined by calcu- 
lating the current at r — > 00 and matching the physical 
boundary current, see more discussions in subsection A. 
below. 

To this point, our problem of finding quantum me- 
chanic solutions (with arbitrary k, v) for monopole scat- 
tering off flux tube have been all set. With these ana- 
lytical solutions at hand, a few discussions are in order 
below. 



A. Scattering Amplitude 

Now we discuss the boundary condition and determine 
the scattering amplitude. As a scattering problem, we 
expect an incident current described by transverse plane 
wave, say e lkx , in the cylindrical setup. Thus we write 
down the asymptotic wave function as 9 



\£k(r — > 00) = e 



ikx 



7? (30) 



Expanding e lkx — e 4fcrcos< A also in terms of e 1 " 1 ^ and 
comparing the above to the large r limit of fff % J (r) from 
eq.(25), we obtain the normalization constant B as 



Bi. 



_ & m{u/2-d) 



(31) 



with the feature \Bk^\ 2 = 1 independent of k, v values 10 . 

The partial-wave scattering amplitude is determined 
via phase shift as 



,i(2Sk, v -dir) _ ^ 



Jmtj> 



(32) 



This gives the partial-wave scattering cross section, or 
more precisely transverse cross " length" , as 



Sk, v = ^ sin (Sk, v 



dir/2) 



(33) 



The total cross section is a sum of the above over all v. 

Examples of 5k, v and S k , v as functions of k for several 
values of v are plotted in Fig. 3. 

Before closing this subsection,we'd like to point out 
that the phase of coefficient given in (31) is related to 
the choice of e lkx as asymptotic incident state (while its 
unity amplitude is general). Physically an incident par- 
ticle can come in from any direction besides x axis, with 
equal probability, so an average over all possible orienta- 
tion of initial k is called for. This can be achieved by first 
doing calculation using (31) and averaging over the <j> de- 
pendence at the end, and the effect of this procedure is 
simply the entire suppression of interference terms among 
different partial waves. 



9 Here we temporarily normalize the incoming current as just 
v = hk/M while in later section additional factor from density 
n will be included. 

10 One should keep open mind in that different boundary 
conditions lead to different weights Bk, v among partial waves. 
It is not clear if there could be choices other than the ones used 
here which can best describe the thermal monopole scattering 
by flux tube field. An extreme example is superconductor 
which picks only B k ,o with all others vanishing. 
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FIG. 3. (color online) (upper) Scattering phase shift 5k,v 
and (lower) scattering cross section Sh,v as a function of kR 
for v — 0(red), v — ±2(blue solid/dashing), and v = ±5(green 
solid/dashing). 



B. Magnetic Current 

The magnetic current generated by single monopole 
during scattering process can be calculated by 



3m 



igh 



2M 



Mc 



(***)C (34) 



J IK 

The nontrivial part is the <j> component: 11 

oo 

J M (k\r)= E iUkMr) (35) 

\fk,u\ 



v— — oo 



gh_ y> 

Mr ^ 



K R 2 ' [ R l 



with 9[x] the unit step function. We may further combine 
±|z/| terms and rewrite it as 



gfi 



j^£&(Mr) (36) 



11 f component is zero and z component is totally irrelevant 
and can also be set to zero by replacing e ** with real sin, cos 
functions. 



d-(l-y-6[l-^}\f k , \ 2 



&(M|r) 

f M {kAr) = v-[\fkA 2 -\h,-»?] 

_d.(l_^).(9[l-I]. [\f k , v \* + \f k ,- v \ a ] 

This expression implies two important points: first, 
significant contribution to magnetic current comes from 
small r part as is evident from 1/r dependence, so partial 
waves with large amplitude at small r (namely small \v\ 
channels) are important; second, according to jff oc v at 
r > R, partial waves with m symmetric to d, namely a 
pair of ±|z/| channels, tend to produce opposite currents 
which substantially cancel each other. It is worth empha- 
sizing that only v = partial wave (the one picked by the 
whole condensate in ANO vertex case, see e.g. [21]) will 
benefit from the first point and at the same time NOT 
suffer from the second point. 

Clearly for each given k the total current jfj(r) to be 
integrated in eq.(18) should be built up from summing 
currents of all partial waves , namely summing j$f(k, v\r) 
over quantum numbers v. One is naturally concerned 
with the convergence of such infinite summation, which 
is basically determined by the large \v\ behavior. We can 
expect that large \v\ partial waves contribute very lit- 
tle to the total current, which bears two simple physical 
arguments: from energy point of view, states with \v\ ex- 
perience centrifugal potential V(r) ~ h 2 \v\ 2 /2Mr 2 while 
the kinetic energy being = h 2 k 2 /2M, so if kR < 
then Ek < V(r — R) which means it is very hard for the 
particle to "climb" up the potential barrier all the way 
into the tube; from the impact parameter perspective, 
states with \ v\ and k have semiclassical impact parameter 
b ~ \v\/k, so if kR < \v\ then b > R which means the inci- 
dent particle will be largely missing the central part and 
thus very little scattered, leading to negligible induced 
currents. This conclusion has been confirmed by exten- 
sive numerical calculation and practically for given kR 
all partial waves with v > l.bkR are vanishingly small, 
as is evident from Fig. 4 to be explained in next subsec- 
tion. 



C. The Total Current 

Now we perform the radial integration needed in 
eq.(18): 

/ 4 ld r= 9 -I(kR) = g -Y,MkR) (37) 

r R 

l v {kR)= / dr&(Mr)/r 
Jo 

In Fig. 4 we plot T v versus kR for various v. The inter- 
esting observation is that the integrated current is nega- 



tive for v < 2, positive for v > 2, while for v = 2 partially 
positive (at small fci?) and negative (at large fci?). This 
result from quantum mechanics perfectly agrees with our 
conclusion from classical treatment in the Introduction 
part, not only qualitatively but even quantitatively: the 
critical angular momentum for current inversion observed 
here L c = u c h with v c = 2 coincides with that predicted 
by eq.(5) once our flux g<&E — d2Tthc with d = 2 is 
plugged in. 




FIG. 4. (color online) Integrated current I„ as a func- 
tion of kR for different values of v: dashed lines are for 
v — O(red), v — l(green), v = 2(blue), while solid lines are 
for v = 3(black), v — 4(red), v — 5(green), v — 10(blue), 
v = 20(magenta), and v = 30(purple). 

Now we perform the final step: namely summing X v 
over v to obtain the integrated total current X. This is 
done numerically, with summation cut v < v cut applied, 
see Fig. 5 for results for various v cu t- As can be seen, 
for the displayed regime kR < 20 the summation is con- 
verged enough as soon as v cut > 20, as the curves with 
v C ut — 20 and v cut — 30 coincide on top of each other and 
are hardly distinguishable. It is this numerically evalu- 
ated function X(kR)(with our highest cut v rut = 30) that 
will be used in subsequent sections. 

The behavior of this function X(kR) has rather non- 
trivial wiggle structure: the general trend is oscillatory, 
with a modest negative part at small kR < 1.42 (basi- 
cally from negative contribution from v = 0, 1) followed 
by a rather high positive peak (dominantly from v = 2) 
between 1.42 — * 2.24. These first two structures, first 
negative then positive, basically cover the interesting re- 
gion of kR(see discussion in next paragraph). Suppose 
there is a flux tube with certain R, then at low temper- 
ature the typical k is small and kR falls within negative 
region which supports the flux tube, while at high tem- 



perature the larger k brings kR beyond the negative re- 
gion into the tremendous positive region which will kill 
the flux tube. So there is a transition with the border 
at kR = 1.42: beyond this point higher partial waves 
with v > 1.5 x 1.42 = 2.13 (which is also close to the 
classical critical value v c = 2) will become dominant. By 
Comparison of this curve with the red dashed one (only 
v = 0), which corresponds to what superconductor can 
do, one understands why a condensate does much better 
in confining a flux tube than a normal thermal ensemble 
can do. The first negative peak, actually the best point 12 
for flux tube formation, locates at 



k m R = 1.076 with X{k m R) = -0.140 



(38) 
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FIG. 5. (color online) Integrated total current I as a func- 
tion of kR for different values of summation cut v cu t'- dashed 
lines are for u cut = 0(red), v cut = l(green), v cut = 2(blue), 
while solid lines are for v cu t = 5(red), v C ut = 10(green), 
Vcut = 20(blue), and u cut = 30(black). 

Finally let's discuss interesting range of kR. Remem- 
ber ultimately we'd like to discuss the flux tube inside an 
ensemble of monopoles with temperature T . So first, the 
tube radius shouldn't be much larger than fic/^r, oth- 
erwise the tube's transverse vibrational modes (w ~ 1/-R) 
get too easily thermally excited, making it unstable. Sec- 
ond, large k should be suppressed by thermal distribu- 
tion, and typical k should be few times fc,eT /Tic. Thus it 



12 One might argue that there will be an even larger negative 
peak at kR ~ 2.77, however to reach that point one requires 
much larger k which usually means broader distribution over 
k around k, and that will easily make the total contribution 
rather small after cancellation with the adjacent large positive 
peak. 



follows that typical values of kR should be of the order 
unity. The evaluated I(kR) up to kR = 20 here should 
be sufficient for later application. 



D. Partial Wave with v — d and Possible Resonance 



general phenomenon and the induced current produced 
by these resonance states actually will spoil the origi- 
nal flux tube field as monopole in such state stays in the 
center of tube and "pushes" field outward rather than in- 
ward: remember the large positive peak in I(kR) (black 
curve in Fig. 5) is precisely due to the contribution from 
v = d partial waves. 





FIG. 6. (color online) (Upper) Phase shift and (lower) scat- 
tering cross section as functions of kR for v = d = 4 which 
show resonance structure, see text. 

The effective potential V e ff with v = d or equiva- 
lently m = is special in that it vanishes at the cen- 
ter r = 0, while for all other v =^ d states there will 
be diverging term ~ (v — d) 2 h /2Mr 2 . The shape of 
it (see Fig. 2) actually indicates possibility for resonance 
to occur. Whether there could be resonance solution 
or not depends on the competition of the localization 
energy and the potential barrier whose peak value is 
E p . = d 2 h 2 /2MR 2 atr = R. To settle this one can look 
at the condition for the wave function (24) to be zero 
right upon r = R (which is very close to the resonance 
situation and gives estimate of kinetic energy). This 
yields the series of particular values of k: k R = 2.576, 
kiR = 5.632, k 2 R = 8.729, k 3 R = 11.847, ... Thus 
clearly to have one resonance level, one needs at least 
E p > h 2 kQ/2M namely d > 2.576. Indeed by fine search 
for resonance structure in scattering phase shift we iden- 
tified one resonance in the case of d — 4, see Fig. 6, with 
k res , m 2.77 very close to the above fco and narrow width 
Tfc ~ 0.1/ R. Nothing similar was found in d = 2. With 
large enough d the occurrence of resonance should be a 



E. Quantum Mechanical Motion of Single 
Relativistic Monopole 

In this subsection we generalize the obtained solutions 
to relativistic case. Now one has to solve Klein-Gordon 
equation (since monopoles are scalar particles) instead of 
Schrodinger equation: 



M 2 c 4 



P 



-C 



* = 



(39) 



Fortunately it turns out that by again writing eigenstate 
of energy e as $ = fk,v(r)e lm ^e %K * z one recovers exactly 
the same radial equation as eq.(22) except for changing 

k = ^(2Me)/h 2 - K 2 to the following 



k = y/(e 2 -AP c 4 ) / '(he) 2 - K 2 



(40) 



So all the exact wave functions obtained in non- 
relativistic case are still solutions to the Klein-Gordon 
equation after the above replacement of k in eq.(24)(25). 
This change of k should be done for all the relevant for- 
mulae above. 

Another important change is for the current equation 
(35): due to relativistic effect the mass M should be 
replaced by e/c 2 , namely 



(e/c 2 



d-(l 



R 2 



).6[1 



R l 



\fk,u\ 2 

(41) 



The same replacement should also be applied to integra- 
tion over current in eq.(37). 

All other aspects remain pretty much the same as in 
non-relativistic case and we skip further discussion. 



V. SELF-CONSISTENT ELECTRIC FLUX TUBE 
SOLUTION 

In this section we will self-consistently determine the 
size R of the electric flux tube carrying flux $e = 
d x (2t{Tic) I g in an ensemble of monopoles with tempera- 
ture T and density n, by using the integrated 13 magnetic 
current obtained in previous section. 



' Since we are not interested in details of the flux tube shape, 
we refrain from doing more complicated local matching of the 



At this stage the issue is to average the integrated total 
current T(kR) over proper thermal distribution n{k\T) 
(through which the medium property comes into play) 
satisfying n = J °° dkn(k\T). We have 



\ jf * jUr = j£-J~dk n{k\T)l{kR) 



(42) 



Note in relativistic case we have to replace the mass M 
by e(fc)/c 2 and move it inside the integration over k. Be- 
low we deal with non-relativistic gas, relativistic gas, and 
optimally correlated ensemble separately. 



A. Non-Relativistic Gas 

In non-relativistic (NR) gas with Mc 2 /ksT large, the 
kinetics are simplified, yet in principle one still needs to 
take into account the quantum statistics, namely using 
the Bose-Einstein(BE) distributions. Only in the non- 
degenerate limit (with monopole gas being not dense) 
one recovers the Boltzmann limit. So we use the BE 



distribution I j(z~ 



n = s x 



(Mk B T 
V 27r?i 2 



I) normalized to density n by 

3/2 



x Li 



-i "- 



(43) 



In the above s is the degeneracy due to internal degrees 
of freedom, fugacity z — e?-l kBT is related to chemical po- 
tential and valued as < z < 1 in NR case, and Li 3 / 2 [z] 
is the polylogarithm function. We then have the n(k\T) 
given by (after integrating out the z momentum) 



n{k\T)dk 



Mk B T 
2ttH 2 



3/2 



x Lii 



2ydy 
(44) 



with the variable y = Uk/(2MkBT). 

Now by combining eq. (18) (42) (44) we obtain the self- 
consistent equation for flux tube size R: 



ttR 2 Mc R 2 Mk B T 



x U 



R 



q = 



h 2 /irMk B T 



(45) 



with the last term U from integration over x = kR 

Lulze-™ 2 ^ 



U[q] 



„ dX X Lis[z] 



l(x) (46) 



current and the \/ x E, as local form of Maxwell equation 
demands. 



The self-consistent equation can be further rewritten 
in an elegant way: 
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(47) 



with \ L = {Mc 2 lg 2 n) x l 2 and X dB = {h 2 /TrM/t^T) 1 / 2 . 
So for given parameters one uniquely determines the flux 
tube size R from the above equation. 

The NR Boltzmann limit, satisfying scale hierarchy 
1/n 1 / 3 » h/y/Mk B T » Ji/(Mc), can be achieved by 

simply replace — - fj^\ in the integration of U[q] 



by e" 



-TTX 2 /2q 2 



Li 3 [z] 



Mathematically this follows from taking 



the z — > limit (with only linear terms left) of both 
polylogarithm functions. 

The results from solving eq.(47) are plotted in Fig. 7. 
Numerically we didn't sec much difference between z — > 
0(blue curve) and z — > I(red curve) limits. As Xl/^b oc 
(T/n) 1 / 2 , the right end of the horizontal axis corresponds 
to high-density/low-temperature regime while the left 
end represents low-density/high-temperature regime. 
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FIG. 7. (color online) R/XdB versus Xl/XcIb from solution 
of eq.(47). The blue curve is for Boltzmann limit (2 -> 0) 
while the red for z = 0.999, see text. 

The distinguished feature is the existence of critical 
point for Xl/Xcib beyond which there will be no self- 
consistent solution: this occurs at roughly the same value 
for both displayed curves and we obtain the following 
condition for the existence of flux tube 



X L 

XdB 



TrM 2 c 2 k B T 
h 2 g 2 n 



1/2 



< 0.13 



(48) 



Physically the above result is very appealing: it 
demonstrates the mechanism of how a flux tube which 
exists in the medium at low T is eventually gone as the 
medium is heated up; on the other hand, for a medium 
with given T it sets up a lower bound of monopole density 
that is required to support the existence of flux tube. 

Another feature is that for each given Xl/ XdB smaller 
than the critical value, there are actually two solutions, 



one with small R/XdB (typically smaller than 1, see 
solid curves) and the other with large R/XdB (typically 
greater than 1, see dashed curves). This is understand- 
able according to the complicated wiggle structure of 
T(kR). The solution with smaller radius is the stable 
one: it is much thinner and thus has stronger electric field 
(E <~ $s/r 2 ), which reflects monopoles more sharply 
near the boundary. The other solution with larger R/XdB 
is unstable and should be discarded. 



B. Relativistic Gas 

In relativistic gas the important scale is set by tem- 
perature, so let's introduce the following dimensionless 
variables: 



R 



w 



hc/k B T 



u = 



Mc 2 
kB~f 



(49) 



The fugacity z, now in range < z < e u , is related to 
density by 



n= s x 



k B T 



3 /-°° t 2 dt/(2w 2 ) 
he J J z -i e V^+t^ _7 

(k B T\ 3 



(50) 



with the number s the degeneracy due to internal degrees 
of freedom. C serves as normalization constant to mo- 
mentum distribution (after scaling momenta by fc^T/c). 
The distribution over k is given by 



n{k\T)dkdK z = s 



kdk dK z /(4n 2 ) 

z -l e e{k,K z )/k B T _ I 



(51) 



with e(k, K z ) = y M 2 c 4 + h 2 k 2 c 2 + h 2 K 2 c 2 . Similarly 
combining the above with eq.(18)(42) one obtains the 
relativistic version of the self-consistent equation 
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We can further rewrite the self-consistent equation as 

2 

(54) 



87r 2 dx \ =w 2 U[w\u,z] 

\ XdB J 



with the newly introduced relativistic parameters Xl 
(kBT/g 2 ^ 1 / 2 and X dB = hc/k B T. 



For given sets of parameters M, n, T (or equivalently 
u, z, T) one can easily find the flux tube size R from 
the above equations by direct numerics. The situation 
is quite similar to the non-relativistic gas which we skip 
further discussion. 



C. Optimally Correlated Ensemble 

Finally let's discuss ensemble beyond an ideal gas. 
Clearly with significant interparticle correlations the en- 
semble may even not be easily describable by any dis- 
tribution, however a typical momentum kr can still be 
invoked. A special situation which we call optimally cor- 
related ensemble is that monopoles from such ensemble 
are largely carrying momenta within very narrow region 
around kx- On the contrary if the ensemble particles' 
momenta are very diffusive in momentum space, it can 
hardly support flux tube. 

In the optimally correlated ensemble, we approximate 
eq.(42) as (assuming NR formulae) 



and the self-consistent equation is then given by 

$ £ gnh - 
= ~ttR 2 = Mc X [ - I{kTR)] 



(56) 



We limit the value of kxR within — 1.42 beyond which 
there won't be flux tube solution, as discussed in Sec- 
tion.IV C. 

The above can be re-organized into 

2d x (k T ■ X L f = {k T Rf x [-l{k T R)\ (57) 

The best situation occurs (roughly) around the negative 
peak in l(kR) given by (38). From this we set a bound 
similar to eq.(48) 



,'Mc 2 k 2 \ 1/2 
k T -X L =[ — 5-21 ) < 0.20 
g z n 



(58) 



VI. DISAPPEARANCE OF FLUX TUBES IN 
QUARK-GLUON PLASMA 

Results from previous sections are general in nature 
and applicable to a variety of plasma physics problems. 
The present section, on the other hand, is dedicated to 
our main application, the physics of sQGP. From now 
on we switch to natural units and systematically put 
h,c,ks = 1. 

The existence of string/flux tubes in the QCD con- 
fined phase T < T c is rather thoroughly investigated on 
lattice, via measurements of static heavy quark poten- 
tials. Static free energy potentials F(T, r) as a function 



of r are only studied for N c — 3 but for number of quark 
flavors Nf = 0, 2 as well as physical QCD, see [32]. Those 
can be used to extract the entropy and potential energy 
separately: the peaks of these quantities (see e.g. Fig. 2 
of [33] ) happen to be exactly at T = T c and then decrease 
toward larger T. The presence of the quasi-linear part of 
the energy and entropy at intermediate r leads to a con- 
clusion that flux tubes still exist at T > T c . Fig. 8 from 
[33] shows how both the internal energy and entropy look 
like at T = 1.3T C . Unlike in the free energy (open squares 
in the upper plot, in which cancellation takes place), the 
internal energy (closed circles) still show at intermediate 
r = (.3 — .7) fm a part linearly dependent on r, while at 
T > 1.3T C it very quickly disappears. 
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FIG. 8. (a) The singlet internal energy, Ui(r,T) (filled 
circles), calculated from renormalized singlet free energy, 
Fi(r,T) (open squares), at fixed T ~ 1.3T C in 2-flavor lattice 
QCD compared to V(r,T — 0) (line) . (b) The corresponding 
color singlet quark anti-quark entropy, TSi(r,T ~ 1.3T C ), as 
function of distance calculated from renormalized free ener- 
gies. 

Why are flux tubes disappearing at large T? It can in 
principle be due to two different changes in QGP happen- 
ing as T grows above T c , to be called (i) electric screening 
and (ii) magnetic penetration. The simplest mechanism 



(i) is that as T grows beyond T c , the density of elec- 
trically charged quasiparticles - gluons and quarks - is 
growing and eventually it becomes large enough to screen 
heavy quarks. The reason for this density growth is the 
decrease in effective masses of electric excitations, which 
are lattice observables by themselves 14 . At very large 
T » T c , in weak (electric) coupling domain, the screen- 
ing of the potential is expected to be described by the 
Debye theory. However Debye theory does not describe 
entropy and internal energy associated with static quarks 
at T = (1 — I.3)T C , even at large distances, as can be seen 
e.g. from calculations of Antonov et al [35]. 

Another effect (ii), discussed for the first time in this 
work, is the penetration of magnetically charged quasi- 
particles (MQPs) inside the flux tubes, which destroys 
them. Indeed, the key parameter kR increases with T 
and reaches the critical value eq.(58) for whether flux 
tube can exist or not. This imposes the following condi- 
tion 



(59) 



Changing T from T c upward the monopoles gets heav- 
ier and their dimensionless magnetic density n/T 3 keeps 
decreasing: eventually this will violate the flux tube con- 
dition. We thus identify the equality in (59) with the 
temperature T m 1.3T C at which local dissolution of the 
flux tubes takes place. 

Furthermore, at T w 1.3T C we expect g 2 /^ir w 1 
[1]. An independent consideration fixes conditions for 
monopole Bose condensation [29] which demands that 
around T c the monopole mass over temperature M/T « 
1-1.2. 

Combining these estimates with our critical condition 
for tube dissolution we obtain the density of magnetic 
quasiparticles at 1.3T C to be 



n MQPs 



(£) 



(60) 



T=1.3T C 



which is within umqps = (4-4 — 6.6) fm~ 3 in absolute 
units. 

Can the density of magnetic objects really be of that 
magnitude (which superficially looks rather high)? This 
estimated density includes in principle contributions from 
all types of magnetically charged objects in sQGP, i.e. 
not only pure adjoint monopoles but also self-dual dyons 
and also dyons containing quarks 15 . 



14 See a related discussion of various color-electric objects' 
effective masses in [34]. There we showed the masses are 
still rather large and their densities rather small at 1-1.5T C , 
limiting the screening. 

5 We recall that monopoles have fermionic zero modes and 
states made of fermions travelling on top of a monopole have 



Let us compare the numbers with whatever is men- 
tioned in literature. We don't know any studies of 
fermionic objects mentioned above. 

Ilgenfritz et al [15] determined their dyon estimate by 
the caloron density, which is reliably calculated from the 
topological susceptibility. After multiplying by 3/2 their 
result for SU(2) we obtain density of self-dual dyons 
to be ridyons ~ 3/m~ 3 . Chernodub and Zakharov [36] 
mentioned the monopole density which is directly es- 
timated from lattice configurations by following gauge- 
fixed monopoles along their trajectories. Their estimate 
is about n mono w 3.5/to~ 3 . The sum of the two is consis- 
tent with the upper end of our estimate of what is needed 
for formation/dissolution of the flux tube. 

Independent comparison can also be made with the 
vacuum (T = 0) monopole density. Bali [22] has mea- 
sured London penetration length by fitting lattice result 
with Abelian Higgs model. From that one can infer the 
monopole density to be as large as 10/to -3 . Bornyakov 
et al [37] gave the vacuum monopole density to be about 
7.5/m~ 3 . All these results are well above our estimates 
for the density at T = 1.3T C "dissolution point". 

We believe all these numbers are consistent and sug- 
gest a coherent picture, of very dense monopole con- 
densate in vacuum, tightly confining electric flux into 
very narrow tubes. When heated slightly above T c the 
monopole condensate changes into a non-condensed en- 
semble of monopoles, which is roughly twice less dense. 
Yet it is still capable of supporting flux tubes survived 
from vacuum, and only around T = 1.3T C the density 
of monopoles drops so low that there won't be flux tube 
any more. At higher T the electric sector becomes more 
and more dominant till eventually small number of heavy 
monopoles become embedded in the perturbative electric 
plasma. 

VII. SUMMARY AND OUTLOOK 

In this paper we have studied stability of the electric 
flux tube in a monopole plasma. Quantum scattering 
of a single monopole on electric flux tube is analyzed in 
great details. Already classical analysis hints on the ex- 
istence of a critical angular momentum dividing the scat- 
tered magnetic currents which support /dissolve the flux 
tube. This finding is quantitatively confirmed by quan- 
tum mechanic calculation, in which we have found exact 
scattering solutions to Schrodinger/Klein-Gordon equa- 
tion in non-relativisitic/relativistic situations. These so- 
lutions allowed us to calculate the magnetic current pro- 



duced, which is then averaged over the monopole ensem- 
ble and used in self-consistent determination of the flux 
tube size. The exact critical condition has been estab- 
lished, and applied to electric flux tube dissolution in 
sQGP system which interests us most. This leads to 
an estimate of total density of magnetic quasiparticles 
nMQPs ~ 4.4 ~ 6.6/m~ 3 at T w 1.3T C , where lattice 
potentials indicate flux tube dissolution. These numbers 
are consistent with other studies using alternative ways 
to estimate magnetic density. 

As mentioned in the introduction, this work is partly 
methodical in nature, ignoring electric quasiparticles 
which would lead to screening and termination of flux 
tubes. The next step we plan to do is obviously inclu- 
sion of both components and calculation of the static 
potentials. Hopefully, when one would consider an ap- 
propriate mixture of electric and magnetic quasiparticles, 
the lattice data on static potentials between electric and 
magnetic 16 charges would be explained. 

In principle, one should go beyond that and calculate 
field distributions around static charges as well. Lat- 
tice studies can be extended to measure directly elec- 
tric/magnetic fields at T <~ T c : in fact the field profiles 
have been measured for flux tubes in vacuum before (see 
e.g. [22]). 

Let us end with the following intriguing question. We 
focused above on electric flux in magnetic media, ignor- 
ing electric quasiparticles and possible dual phenomenon 
- a magnetic flux tube in an electric plasma. (We only 
mentioned their existence at low T high density regime, 
in a color superconductor.) Now, may somewhere along 
the electric-magnetic equilibrium line there be conditions 
supporting stable flux tubes of both types at the same 
time? It is known that confinement of both is impos- 
sible, but in a uncondensed plasma regime it may still 
be the case. A natural place to look for a QGP with 
intertwined electric and magnetic flux tubes is at T less 
or of the order of T c , close to the place where three ma- 
jor phases - hadronic, color superconductor and QGP - 
meet. Although it is quite challenging task to get into 
this region using lattice gauge methods, the task is not 
hopeless. 
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to be included as well. In supersymmetric theories those form 
spin-1/2 and even 1 magnetic objects, which are needed by 
supersymmetry to produced appropriate supermultiplets in- 
cluding the usual scalar monopoles. 



Those are given by the expectation value of the so called 
't Hooft loop: we have not discussed them in this work. 
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